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Simulations of supersymmetric Yang-Mills theory 
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Results of a numerical simulation concerning the low-lying spectrum of four-dimensional JV = 1 
SU(2) Supersymmetric Yang-Mills (SYM) theory on the lattice with light dynamical gluinos are 
reported. We use the tree-level Symanzik improved gauge action and Wilson fermions with stout 
smearing of the gauge links in the Wilson-Dirac operator. The configurations are produced with 
the Two-Step Polynomial Hybrid Monte Carlo (TS-PHMC) algorithm. We performed simulations 
on lattices up to a size of 24 3 • 48 at /3 = 1.6. Using QCD units with the Sommer scale being 
set to ro = 0.5 fm, the lattice spacing is about a ~ 0.09 fm, and the spatial extent of the lattice 
corresponds to 2.1 fm to control finite size effects. At the lightest simulated gluino mass our 
results indicate a mass for the lightest gluino-glue bound state, which is considerably heavier than 
the values obtained for its possible superpartners. Whether supermultiplets are formed remains to 
be studied in upcoming simulations. 
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1. Introduction 

Supersymmetric theories have received increasing interest in elementary particle physics. The 
supersymmetric extension of the Standard Model with <jf = 1 supercharge and other models be- 
yond the Standard Model have supersymmetry (SUSY) as an essential ingredient. 

The minimal supersymmetric extension of the S\J(N C ) gauge theory describing self-interactions 
of gauge fields A^, corresponding to the gluons (g), is given by the jY = 1 Supersymmetric Yang- 
Mills (SYM) theory. The supersymmetric partners of the gluons are described by spin- 1/2 Ma- 
jorana fermion fields X a {a = 1 . . .N% — 1), the gluinos (g). Compatibility of SUSY with gauge 
invariance requires that the gluinos transform in the adjoint representation of the gauge group. This 
theory describes the interactions between gluons and gluinos. The Lagrangian of SYM theory in 
the continuum, including a SUSY breaking mass term, reads 

! i _ mi - _ 

^sym = -^v^ v + ^"Wf ( v wr - -f xar ■ (L1) 

The mass term here introduces a soft breaking of supersymmetry. 

In the low-energy regime, where the interactions become strong, arguments based on the low- 
energy effective Lagrangian approach [|l|, ^ predict the occurrence of non-perturbative dynamics 
like confinement and spontaneous chiral symmetry breaking in SUSY gauge theories. Confinement 
is realised by colourless bound states. In the case where the last term in Eq. ( |iH| ) is switched 
off {nig = 0), the anomalous global chiral symmetry U(l)^ is present. The anomaly does not 
break the global chiral symmetry completely and a discrete subgroup Z2n c remains. As in the case 
of QCD, the discrete chiral symmetry is expected to be spontaneously broken to Z2 by the non- 
vanishing value of the gluino condensate (A A). The consequence of this spontaneous breaking is 
the existence of N c degenerate ground states with different orientations of the gluino condensate. 

Another interesting aspect of SYM is its equivalence to QCD with a single quark flavour 
(Nf = 1 QCD) in the limit of a large number of colours (N c — ► 00), where the Majorana spinor 
is replaced by the single Dirac spinor The latter model is also object of investigation by our 
collaboration [H]. 

A lattice formulation of SYM suitable for numerical simulations, employing Wilson fermions, 
has been proposed by Curci and Veneziano [g]. First non-perturbative investigations of SYM on 
the lattice using this formulation have been performed by the DESY-Munster-Roma collaboration; 
for a review see Ref. |@] and references [|, g, g, |^ 01. SUSY is broken explicitly by the lattice 
discretisation. In the Wilson approach the mass term and the Wilson-term break both chirality and 
SUSY. Both symmetries are expected to be recovered in the continuum limit by tuning the relevant 
bare mass term to its critical value corresponding to a massless gluino (trig = 0), and the gauge 
coupling towards zero. 

In recent years, simulations of jY = 1 SYM on the lattice using Ginsparg-Wilson fermions 



with good chiral properties, such as domain wall fermions, have been initiated Q12| , |13j, [14p. For 
large lattice volumes and small lattice spacings these formulations require, however, a significantly 
larger amount of computing resources than the Wilson formulation. The gain of no need for tuning 
the position of the zero gluino mass point does not compensate by far the advantage of Wilson 
fermions. 
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Here we report on recent work continuing the project of the DESY-Munster collaboration for 
the simulation of J/ = 1 SU(2) SYM. We present results in the light of the newly used TS-PHMC 
algorithm and improved actions. 

The most important characteristics of the theory is the mass spectrum of bound states, for 
which the low-energy effective theories predict a reorganisation of the masses in two massive Wess- 
Zumino supermultiplets at the SUSY point [jl], ^J, where the soft breaking vanishes. The introduc- 
tion of a small gluino mass removes the mass degeneracy between the supermultiplet members. 



2. Lattice formulation of jY = 1 SYM theory 

For the gauge fields we employ the tree-level improved Symanzik (tlSym) action. The gluinos 
are represented by Majorana fermions X a in the adjoint representation. The fermion part of the 
Curci-Veneziano action, describing the gluinos, is given by 

^^^I^W^W-fll^^+A^WCi+rM^W+^W^WCi-rM^^+A)], (2.1) 

^ x x H 

where K is the bare hopping parameter which encodes the bare gluino mass K = (2mg o + 8) ~ 1 ■ The 
real orthogonal matrices V^(x) are the gauge links in the adjoint representation. 



The links U xfl in the fermion action can be replaced by sfowf-smeared links [15]. This has 
the advantage that short range topological defects of the gauge field and the corresponding small 
eigenvalues of the fermion matrix are removed. We prefer to keep the action well localised and 
hence only perform a single stout-smearing step. 

Similarly to QCD, the mass term proportional to mgo breaks chirality. In the present case it 
also breaks the supersymmetry. A massless gluino, trig = 0, is obtained by tuning the relevant bare 
mass term to its critical value (mo — > mo cr ) or equivalently K — > Kc r . 

The fermion action can be rewritten in term of an antisymmetric matrix M = ^€Q, where Q is 
the non-hermitian fermion matrix or lattice Wilson-Dirac operator for Dirac fermions in the adjoint 
representation, and *<£ is the charge conjugation matrix in the spinorial representation. Integration 
of the fermionic variables yields the Pfaffian of M, whose absolute value equals the square root of 
the fermion determinant. Effectively, this corresponds to a flavour number Nf = I. In the Wilson 
setup the Pfaffian can become negative even for positive gluino masses. 

In our numerical simulations we include the dynamics of the gluino by the Two-Step Polyno- 
mial Hybrid Monte Carlo (TS-PHMC) [16] algorithm with flavour number Nf = This has the 



consequence that only the absolute value of the Pfaffian is taken into account in the updating of the 
gauge field configuration. The sign of the Pfaffian has to be included in a reweighting step when 
calculating expectation values. It can be shown that the sign of the Pfaffian is equal to the sign of 
the product of half of the doubly degenerate negative real eigenvalues of Q. For positive gluino 
masses sufficiently far away from zero a negative sign of the Pfaffian rarely occurs in the updating 
sequence and therefore a sign problem does not show up. 

The TS-PHMC algorithm turned out to be very efficient in producing short autocorrelations 
among the gauge configurations. For instance, in the stout-smeared runs on a 24 3 • 48 lattice the 
integrated autocorrelation of the average plaquette (which belongs to the worst quantities from the 
point of view of autocorrelations) did always satisfy %f^ q < 10. 
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The values of the gauge coupling parameter j8 can be fixed by investigating the static potential 



of an external fundamental charge and extracting the Sommer scale parameter ro/a Jl7[] . Note that 
in analogy with QCD, we set the value of ro by definition to ro = 0.5 fm. In this way we can use 
familiar QCD units for physical dimensionful quantities. 

The TS-PHMC runs were done on 16 3 • 32 and 24 3 -48 lattices at j8 = 1.6 at various values of 
the hopping parameter K. On the 24 3 • 48 lattice four points have been simulated with stout-links. 
The lattice spacing amounts to a ~ 0.09fm. The lattice extension L ~ 2. 1 fm is expected to be large 
enough to allow control over finite volume effects on the bound state masses. 

The dimensionless quantity M r = (rom^) 2 , where m n is the pion mass in the corresponding 
theory with two Dirac fermions in the adjoint representation, is expected to be proportional to the 
gluino mass (see below), and can be considered to be a measure of it. 

An issue in lattice simulation is the lightness of the dynamical quarks which leads to slowing 
down of the update algorithms. With the new TS-PHMC algorithm the lightest adjoint pion mass 
in our simulations was about 440 MeV. Simulations for smaller gluino masses and/or finer lattice 
spacings are going on presently. 

3. Low-lying bound state spectrum 

For the investigation of the low-lying bound state spectrum we concentrate on the projecting 
operators employed for the construction of the low-energy Lagrangians of [|I|] and These are 
expected to dominate the dynamics of SYM at low energies. Previous experience on the field can 
be found in [f7|] and [g]. We investigate spin-0 gluino-gluino bilinear operators (adjoint mesons), a 
spin- 1/2 mixed gluino-glue operator and spin-0 glueball operators. In some cases smearing tech- 



niques such as APE Q18| ] and Jacobi smearing [ |19| ] are applied in order to increase the overlap of 
the lattice operator with the low-lying bound state. 

The adjoint mesons are colourless composite states of two gluinos with spin-parity + and 0~ . 
In analogy to flavour singlet states in QCD we denote the former a-r\' and the latter a-fo, where the 
prefix a stays for "adjoint". The associated projecting operators are the gluino bilinear operators 
&mes = where T = 75, 1 respectively. 

The resulting propagator consists of connected and disconnected contributions. The exponen- 
tial decay of the connected term defines the adjoint pion mass m a - n . This quantity, even if not 
associated to a physical state of SYM, can be used to characterise the gluino mass. Indeed, accord- 
ing to arguments involving the OZI-approximation of SYM, the adjoint pion mass is expected to 
vanish for a massless gluino and the behaviour m 2 . n oc nig can be assumed for light gluinos [jl], g]. 

For the positive parity glueball + we adopted the simplest interpolating operator built from 
space-like plaquettes. In order to improve the signal we applied APE smearing with the variational 
method [§(J. 

Also here, as for the scalar a-fo, the gauge samples generated with stout links turn out to 
generally give better results for the glueball masses. 

The gluino-glueballs (g-g) are spin-^ colour singlet states of a gluon and a gluino. They are 
supposed to complete the Wess-Zumino supermultiplet of the adjoint mesons We adopt for this 
state the lattice version of the gluino-glue operator Tr c [FaA] |JTJ|, where the field-strength tensor 
Ffiv(x) is replaced by the clover-plaquette operator P^, v (x) pl|,|9|]. We apply APE smearing for the 
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links and Jacobi smearing for the fermion fields in order to optimise the signal-to-noise ratio and 
to obtain an earlier plateau in the effective mass. 

A vanishing gluino mass is a prerequisite for supersymmetry in the continuum limit. Therefore 
the point corresponding to a massless gluino is of high interest in lattice simulations of SYM. With 
Wilson fermions this point must be located by tuning the hopping parameter. The subtracted gluino 
mass can be determined in a direct way from the study of lattice SUSY Ward-Identities (Wis) [21] 
and, in an indirect way, from the vanishing of the adjoint pion mass. Indeed, as mentioned above, 
the pion mass squared (a«%) 2 is expected to vanish linearly with the (renormalised) gluino mass. 
Both the Wis and adjoint pion mass methods give consistent estimates of the critical hopping 
parameter K cr corresponding to vanishing gluino mass. 

The masses of the lightest bound states of low-energy J{ = 1 SYM determined in this work 
are shown in Fig. [|. The masses of a-fo and the glueball from the unstout ensembles have a very 
low signal-to-noise ratio and are not shown. The results at the smallest gluino mass are preliminary. 
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Figure 1: Low-lying bound state masses of JV = 1 SU(2) SYM as a function of the 
adjoint pion mass squared (physical units). 

The masses in Fig. [j] are multiplied by the Sommer scale parameter and plotted as a function 
of the squared adjoint pion mass. The lightest simulated adjoint pion mass is about 440 MeV in 
our units. The vertical line in Fig. |] highlights the massless gluino limit where SUSY restoration 
is expected up to &{a) effects. 

The bound state masses appear to be characterised by a linear dependence on {rom n ) 2 . The 
gluino-glueball (gg) with a mass of about 1440 MeV (in our unit where ro = 0.5fm) turns out to 
be considerably heavier than the a-r]' with a mass of 810 MeV. Furthermore, the masses of the 
scalar glueball and the scalar meson a-fo converge to a common point with the pseudoscalar near 
the region where SUSY is expected. The behaviour of scalars is compatible with mixing between 
+ glueball and a-fo. The behaviour of masses suggests a lower supermultiplet, while the spin- 1/2 
candidate remains heavier up to the smallest simulated gluino mass in this simulation. Whether this 
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outcome is a discretisation artefact or a physical effect, as claimed in [22], should become clear in 
future studies at finer lattice spacings. As the data at small gluino mass are preliminary, it would 
be premature to draw conclusions. 



4. Spontaneous symmetry breaking 

An important feature of the strong interaction dynamics of SYM theory is the spontaneous 
breaking of the discrete chiral symmetry. The appearance of two degenerate ground states in case 
of SU(2) gauge group was observed in [ p3| ] on a small (6 3 • 12) lattice by neglecting the sign of the 
Pfaffian. 

Recently we repeated this computation on a larger (16 ) lattice and also taking into account 
the Pfaffian signs. The observed distribution of the gluino condensate is shown in Fig. ^[ 
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Figure 2: Distribution of the gluino condensate at j3 = 1.6, IC= 0.158. 

The distribution cannot be fitted with a single Gaussian, but a two-Gaussian fit shown on the 
figure, corresponding to the two ground states, is reasonable. From the difference of the positions 
of the two peaks one can, in principle, determine the magnitude of the gluino condensate. This 
however requires the knowledge of renormalisation constants and therefore is left for a future work. 




5. Conclusion 

In this work first quantitative results of the low-energy spectrum of JV = 1 supersymmetric 
Yang-Mills theory are obtained. Physical volumes larger than 2 fm have been simulated, which 
is the volume usually required for spectrum studies in lattice gauge theory. The comparison of 
masses on different volumes in otherwise same conditions reveals negligible finite size effects 
at least for moderate gluino masses. We have collected higher statistics and have used efficient 
dynamical algorithms such as TS-PHMC, which is suitable for light fermion masses. In addition, 
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the supersymmetric Ward identities and other observables like the confinement potential and the 
gluino condensate have been investigated. 

From the results of the mass spectrum the question of the gluino-gluino and gluino-glueball 
mass splitting remains open. It can only be answered by further simulations allowing an extrapola- 
tion to the continuum limit. 
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